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1 Introduction 

Recent observations [1] strongly support that primordial inflation is the theoretical expla- 
nation of how the currently observed Universe was formed at the early stages. Alongside 
with observations theoretical approaches also show how nice inflation can be connected to 
the nucleosynthesis and subsequent appearance of the particle standard model. A number of 
inflationary scenarios are reviewed in [2] and references therein, for instance. 

Even though inflation is a great model for many reasons it has problems one of which is 
the lack of the UV completion. To be more precise it is not UV-complete in the framework 
of the Einstein's General Relativity (GR) since geodesies are not past-complete. This is a 
general statement and it is known as the Big Bang singularity elaborated in [3-5]. One can 
find that alternatives to Big Bang such as "emergent" Universe or bouncing Universe [6] 
hit the singularity theorem by Hawking and Penrose [7] as long as we are in GR and the 
space-time is of the Friedmann-Lemaitre-Robertson-Walker (FLRW) type. 

One of the possible resolution is a modification of GR. This can be done in general 
in a number of ways and one may have an insight in this using the review paper [8] and 
references therein, for example. It is inevitable that any modification of gravity introduces 
higher derivatives and only special structures like Gauss-Bonnet term or Lovelock terms in 
general [8, 9] preserve the second order of the equations of motion but this is applicable only 
in more than 4 dimensions. On the other hand finite higher derivatives lead to ghosts due 
to the Ostrogradski theorem [10]. Having all orders of higher derivatives may open a way 
to evade the Ostrogradski statement and a successful attempt in this direction was made 
considering a special class of gravity modifications where higher curvature corrections are 
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accompanied with non-local operators [11, 12]. Analysis in those papers shows how one can 
construct a ghost-free and asymptotically free modification of GR featuring a non-singular 
bouncing solution, and which resembles the GR itself in the IR limit. . The further study 
of this model [18] has shown that the model features the expected perturbative spectrum at 
late times and is stable with respect to small isotropic inhomogeneous perturbations during 
the bounce phase and in parallel a more general model [19] (see also [20]) was proposed and 
considered in the Minkowski background. 

Absolutely the non-local operators is what makes these models novel and the operators 
in question are of type of analytic functions of the covariant d'Alembertian operator, i.e. 
^(D). 2 We find it is quite important that the initial inspiration for introducing these opera- 
tors came from string field theory (SFT) models which as the whole theory is a UV-complete 
non-perturbative description of strings. We refer the reader to more stringy oriented lit- 
erature [23-28] for a more comprehensive overview of this aspect. A decent progress was 
achieved in studying non-local scalar field models derived from SFT in the cosmological con- 
text [29]- [32] as well as exploration of other aspects of this type of models including their 
thermodynamics [33-35]. The major question of rigorous derivation of the modified GR ac- 
tion involving the non-local operators of interest from the scratch, i.e. from the closed SFT, 
is still awaiting for its resolution but this is beyond the scope of our present study. 

The present paper is aimed at extending and continuing both papers [18, 19]. In [18] 
a lot of technical issues were solved for a model which contains the scalar curvature squared 
non-local term. The main focus there is perturbations around a bouncing solution. In [19] 
non-local terms containing the Ricci and Riemann tensors squared were added but only 
analyzed around the Minkowski background. We literally want to join those works in the 
present paper. We confine ourselves by considering the FLRW type of the metric and positive 
cosmological term A. Having these in mind we focus on deriving full equations of motion and 
bringing them to a form one can use in the future study. Our main goal is to present here 
classical solutions and develop the perturbation technique at least in some regimes or around 
some backgrounds. This is crucial for claiming the model is viable or not for the purpose of 
describing the non-singular bounce. 

The paper is organized as follows: In Section 2, we formulate the model and introduce 
the relevant notations. Also we write down full equations of motion relevant for the FLRW 
type metric. In Section 3, we demonstrate explicitly that one of the known bouncing solutions, 
namely the cosine hyperbolic, is the solution to our extended equations of motion as well. In 
Section 4, we perform an attempt to construct more solutions and develop relevant techniques. 
In Section 5, we are building the road to analyzing perturbation. We find out we are lucky to 
draw the stability statement about our exact bouncing solution from general arguments as 
well as sketch the derivation of the closed system of perturbation equations for the de Sitter 
asymptotic. In Section 6, we summarize what is done and formulate open questions for the 
future study. 

1 Similar approaches involving non-local models were used in other cosmological and gravity contexts in 
the literature [13-17]. 

2 Theoretically motivated operators such as I/O were considered, for instance, in [21, 22] and references 
therein. 
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2 Action and equations of motion 



We focus on the model described by the following non-local action 

S = J ' d 4 x^{^R+^{RF 1 {a)R + RZT 2 (v)R» + C^ a pT 4 (n)C^ a P)-A\ (2.1) 

where we limit ourselves with 0{R 2 ) corrections. Here the dimensionality is manifest and in 
the sequel all the formulae are written having 4 dimensions in mind, Mp is the Planckian 
mass, A is a cosmological constant and A is a dimensionless parameter measuring the effect 
of the 0(R 2 ) corrections. The most novel and crucial for our analysis ingredients are the 
functions of the covariant d'Alembertian operator J 7 /. For simplicity to avoid extra compli- 
cations we assume that these function are analytic with real coefficients // n in Taylor series 
expansion Ti = J2 n >o fin^ n /M* n - The new mass scale determines the characteristic scale of 
the gravity modification. We assume it universal for all Ti and refer the reader to [11] for a 
detailed discussion of this new physics parameter. Also apart from the canonical usage of the 
Riemann tensor we use the Weyl tensor C£ v p which is coming from the Ricci decmposition 

In this formula we use slightly unusual position of indexes which is useful in performing further 
computations. The reason to use the Weyl tensor is because C£ v p = on a conformally flat 
manifold which is the case for the FLRW metric. We are focused on the FLRW cosmologies 
and thus will benefit out of this. Indeed, it means that the Weyl tensor squared does not 
show up in the background at all and only becomes relevant in perturbations. Moreover, even 
in perturbations the only non-vanishing contribution is the one where both Weyl tensors are 
perturbed and the non-local functions J-4 takes its background form. 

This action appears in [19] 3 and is the most general covariant non-local action up to the 
square in curvature terms and analytic operator functions J-j. 4 Furthermore we note that 
working in 4 dimensions we can assume f 2o = because using that the Gauss-Bonnet scalar 
is a total derivative and combining this with the Ricci decomposition we can have only R 2 
when no d'Alembertian operators are in between. In other words among the terms without 
d'Alembertian operator insertions only R 2 survives on the FLRW backgrounds. This does 
not work for non-constant terms in Ti though. 

Equations of motion for action (2.1) can be derived by a straightforward variation and 
are as follows 

[Mp + 2\T 1 (\3)R]GZ = TP - ASH - ^RTi(a)Rd^ + 2A(VU - 6^)Ti{a)R- 



-2XR^ 2 (a)R^ + ^R^ 2 (a)R^ 



( 11 \ < 2 - 2 ) 

+2A f V p V v T 2 iU)RW _ _ n jr 2 ( D )^ _ -^V CT V p T 2 ip)R ap J + 

-f A/Ci£ - \% {IC^ + £1) + A/C 2 £ - ^ {ICC + IC 2 ) + 2AA£ + 2AC£ 



3 It is eq. (27) there. We have however less terms because possible contractions of the covariant derivatives 
with the Ricci tensor can be eliminated thanks to the Bianchi identity or converted to the higher order in 
curvature terms using the commutation relations for the covariant derivatives. 

4 Also questions of a ghost-free gravity modification are addressed in [36] using different action still using 
non-local operators constructed out of the d'Alembertian operator. 
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where we have defined: 



oo n—1 oo n—1 

n=l 1=0 n=l Z=0 

oo n—1 oo n—1 

^ = E^E v^^v^-^f , £ 2 = E ^ E ^ (0 ?« M t 

n=l «=0 n=l Z=0 

oo n—1 

a* = E/a„E v ^ iJ( ' ) ? VMiJ(n "'" 1) - - v^ B (»-i-i)3] 

71=1 Z = 

and is the part coming from the variation of the Weyl tensor squared piece. As far as 
we mostly concerned about FLRW cosmologies we simply neglect this contribution in the 
background as it is zero. We shall return to it later during the consideration of perturbations. 
Here G# is the Einstein tensor, R^> = D n R, R^ = n n R^, and Ttf is the matter stress 
tensor if one is present in the system. 

One would benefit considering first the trace equation 

—MpR = T — 4A — 6Anj r i(n)i? — \{K,\ + 2K,\) — 
-\UF 2 {u)R - 2XV p V^ 2 {a)R^ J - A(/C 2 + 2JC 2 ) + 2AA + 2AC 

where quantities without indices denote the trace and significant simplifications are obvious. 



3 Exact solutions based on the recursion relations 



3.1 Recursion relations 

In what follows in this Section we use quite intensively ideas and the knowledge accumulated 
in [18] (see also [12]) where this model without T 2 and T± pieces was considered and solutions 
(and their construction) were discussed. We therefore omit extensive citing of these papers 
and do refer the reader to those manuscripts for a detailed analysis of a model which is 
T2 = Fa = in our notations. 

Clearly those equations have as solutions Minkowski (if A = 0) and de Sitter (if A = 
const > 0) backgrounds. The question is however whether we can find other solutions 
especially focusing on the possible bouncing backgrounds. Equations are definitely highly 
non-trivial but a decent progress was achieved in [11] by employing the ansatz 

aR = nR + r 2 (3.1) 

in the absence of J- 2 and JF^. This makes the system much simpler creating the recursive 
relation for and one can try to find exact solutions. Known exact solutions (apart from 
M4 and dSi) include 

a(t) = aocosh(crt) found in [11] (3-2) 

and 

a(t) = aoexp (^ 2 ) found in [37]. (3-3) 
where a(t) is the scale factor of the FLRW metric g^ v = diag(— 1, a(t) 2 , a(t) 2 , a(t) 2 ) . 
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For the Ricci tensor term we suggest the following relation as a quite general ansatz 



UR $ = E (n m Rff + T2 m 9^R (m) ) + rg^ (3.4) 

m=0 

where T\p and r are constants. Surely we require n to be finite. This would manifest the 
recursion relation when acting by extra D-s and give us a chance to have this as an exact 
solution to the equations (2.2). Needless to say that smaller n makes the succeeding analysis 
simpler. 

The explanation why finding of such a recursive relation is almost enough is based on the 
following arguments. One should start with the trace equation (2.3). If there is no matter or 
only radiation (which is conformal and traceless) is present than the trace equation contains 
only the gravitational sector of the model. It can be solved simply by imposing null conditions 
on coefficients in front of different powers of the curvature. These are some conditions on 
the model parameters but clearly the restrictions are very weak as far as n is finite. Further 
we must check one of the component equations (we shall do for (00) component) in (2.2) 
whether it is satisfied or not. But thanks to the Bianchi identity constraint any discrepancy 
we can meet would be of the form of the radiation energy density, i.e. ~ a(t)~ 4 . Hence the 
only question would be how much radiation must we inject in the system and would it have 
good or ghost sign of the energy density. 

Of course the proposed above ansatz is of extreme generality and in order to be more 
specific we take its following simplified version into the consideration 

U n Gi = a n u& v + p n G^ n > 2, where G£ = R» - -^R (3.5) 

The recursion relation for n = 2 can be easily solved 5 to give 



"2 ± + 4/3 2 

a n = c + s + + cs_, s± = y — , p n = a n -ip 2 (3.6) 

and c± are to be found examining explicit values of a^, (5$. However there is a limiting value 
n = 4. n < 4 gives a hypothetical chance to solve the recursion relation while greater n will 
result in general in algebraic equations of a degree higher than 4. 

3.2 Exact analytic bounce 

This choice of the simplified ansatz is motivated by the fact that for solution (3.2) we have 
after some simplifications 

a 2 G£ = Ua 2 nG£ - 40a 4 G£ (3.7) 

Below we mainly study the particular solution (3.2) meaning that we have specific 
s±, c±. One can, surely, track all the steps considering (3.5) in general without specifying 
particular values for a 2 and This is interesting in case more solutions satisfying this 
ansatz can be found. Going general one would come to a bit more complicated expressions 
than we will. We do not do this at the moment for the sake of clarity. Moreover (3.5) is 
already the simplified version of (3.4) which is worth to analyse if a really general scenario 
is of interest. We thus focus mostly on analysing found solution (3.2). 



3 We thank at this point Dan Thompson for illuminating this issue. 
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On solution (3.2) we find 

«™ = ^2 (*1 " Pn = -^fisr 1 ~ »2 _1 ), *1 = !0^, «2 = 4a 2 . 

Note that this relations can be used for n = 1 and n = as they generate correct values 
a\ = 1, Pi = and ao = 0, /?o = 1 respectively which are reasonable. Even though 
/20 = as explained in the previous Section we shall use this property formally to simplify 
the evaluation of the sums in the equations of motion. Apart from this the first ansatz (3.1) 
is also satisfied with r\ = 2a 2 , r 2 = —24a 4 . We thus have 

n n R = r?(R + r 2 /n), □"GJJ = ^ (u&i - s 2 Gt) - (u&i - Sl & u 

These relations yield 

K$ =F£\ri)d»Rd v R, 
ICi = ri F?{ri)R 2 + 2r 2 F?(n)R - {F x {r x ) - fi )R^+ 

+ ijpW (ri) _ (iFl(ri) _ /lo)(r2/ri) 2 ) 
IC 2 =S 1 4 1 \s 1 )S 1 a pS£ + S 2 F { 2 l \ S2 )S 2 p 2 i - • F2(Sl) ~- F2(S2) (s 1 + S^S^ 

S\ — s 2 

+\ (nF { 2 l \n)R 2 + 2r 2 4 1) (r 1 )R - F 2 (n)R^ + ^^(n) - F 2 ( ri )(r 2 / ri ) 2 ) , 
K =V /3 [^ 1) ( Sl )(5 1 ^V^5 1 3 - V^S^SC) + F { 2 \s 2 ){S 2 ^S 2 Z - VSfiS-Z)- 

- • Fa(si) ~ Jc ^ S2 \ s 1 ^s 2 i - VSxfal + S 2 ^StZ - vstfsd)], 

si — s 2 ' ' ' 1 

Fx{U)R = Ji(n)i2 + (Ji(n) - fo)-, = F 2 ( Sl )S^ - F 2 (s 2 )S 2 £ + h»F 2 (a)R 

n 4 

where we have defined 

not - s 2 & v At Dg - Sl G(^ 

liy 6a 2 ' 2l/ 6a 2 

The slightly counter- intuitive remixed position of subscripts 1,2 is correct here. F X2 is the 
derivative with respect to the argument. In order to cancel in the trace equation structures 
containing and (DG^) 2 we have to impose the following conditions 

} ( S i) = J$\s 2 ) = Jb(si) - J- 2 (s 2 ) = (3.8) 
and to cancel structures containing (d^R) 2 and R 2 

^i (1) (n) + ^ 1) (n) = 0. (3.9) 
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All other equations are also simplified on imposing the above conditions. In particular we 
have 

=M*i)(GZ + \KR) + \s^ 2 ( n )(R + r 2 /n) 

or equivalently 

Careful substitution in the trace of the Einstein equations (2.3) reads 



(3.10) 



(3.11) 



A 



-M 2 P R = T - 4A - 6A(Ji(n) + -^{r^nR + r 2 ) - ±T 2 ( Sl )( n R + r 2 )+ 

+2A(7i(n) + - A T 2 {n) - h )-(R + r 2 /n) 

4 r\ 
and we solve it (assuming T = 0) by imposing 

^ - \^{si) - 6A(7i(n) + \Mn)) + 2A(^(n) + ^(n) - / l0 )!J = o 

_Mp_r 2 

4 n 

These relations in terms of a are as follows 

Af 2 \ q 

_£ _ -J- 2 ( 4f x 2 ) - 18A^(2cx 2 ) - -\F 2 {2a 2 ) + 12A/ l0 = 

3MJ.O- 2 = A 

Even though some tuning is required for functions T\ 2 the solution to the above conditions 
is really ambiguous. 

Substituting the obtained results in (00) component of (2.2) and performing a num- 
ber of manipulations we find that only terms proportional to a constant, 1/ cosh (at) 2 and 
1/ cosh(<ri) are present and moreover upon applying the above conditions (3.12) one is left 
with the following combination 

* = (* w + i^ 2(ri) + ^ 2<S2) ) (3 ' 13) 

where we clearly see that we are really left with terms resembling the radiation (i.e. ~ a(t)~ 4 ). 
What is remarkable here that there are two possible resolutions to have the equations solved 
completely. 

First way is to avoid extra radiation requiring 

Mn) + + = o (3.14) 

This means that extra parameter adjustment is important which however does not seem to 
be unnatural. Equations do not trivialize in this case and we can analyse such a configuration 
consistently. Note, it is not the case if J- 2 = (see [18] for details). 
Second way is to admit some amount of radiation with a weaker requirement 

^i(ri) + i-F 3 (n) + ^2{s 2 ) < (3.15) 
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so that this radiation has the positive energy and is not a ghost. 

We therefore see that solution (3.2) equipped with conditions (3.8,3.9) features a number 
of very nice properties: (i) it has non-singular bounce, (ii) it comes at large time to the de 
Sitter phase rather than eternal superinflation, (hi) it requires just few fine-tunings in the 
system which still leave us with an enormous freedom and (iv) it does not have to be supported 
by extra sources if we impose an extra condition (3.14) (to avoid this extra tuning we would 
assume (3.15) to avoid ghosts in the system). 

4 More analytic cosmological solutions? 

In this Section we provide several solutions which unfortunately in most cases do not con- 
tribute to the question of the non-singular bounce but are nevertheless interesting for the 
future developments in this type of theories. 

4.1 Closely related solutions 

Apart from solution (3.2) analyzed in the previous Section one can find after simple analysis 
that the following are also analytic solutions: 



They can be quite simply obtained using the already known one and unfortunately none 
of them contributes to our main subject of a non-singular bounce. Indeed 

• Solution (4.1a) is clearly not a bouncing solution and we mention it here just for the 
sake of completeness. 

• Solution (4.1b) is just our already explored solution (3.2) with an imaginary parameter 
a. It has periodic zeros for the scale factor. It can be used as a test-bed for discussing 
possible stages of the Universe evolution if they obey cos (at) behaviour. 

Obviously there should be other solutions in the system. At first let us mention that 
(3.3) does not generate a recursive relation for the Ricci tensor terms. One could think also 
that there is a possibility to solve equation (3.10) just as a differential equation without 
adjusting the coefficients but this turns out to be incompatible with the ansatzes. Indeed, 
equation (3.10) was obtained using (3.1) and (3.5) while a solution to it does not satisfy these 
ansatzes identically. Even just finding a solution other then presented in the previous Section 
to the ansatzes conditions is not an easy task as it comes to solving at least the third order 
non-linear differential equation. As of now it is an open question how to construct more 
solutions in this model using some ansatz (see [38] for various proposal regarding solving the 
ansatz condition). 

So on one hand we absolutely understand that taming all the non-localities by virtue 
of a recursion relation is the good technical point of the all considered above solutions. This 
also must help in studying perturbations around those particular backgrounds. On the other 
hand however this significantly reduces the possible solutions. We thus want to go beyond the 
recursion relations and to find a more general approach to the problem of solving equations 
of motion. 



a(t) = ao sinh(<7i) 
a(t) = ao cos(at) 



(4.1b) 
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4.2 Model reformulat ion using Gy 

The first useful technical step is a passage from the Ricci tensor to the traceless analog of the 
Einstein tensor Gy. We mention in this regard that appearance of a combination T\ + \T 2 
is not spontaneous because we can rewrite the initial action (2.1) in terms of Gy as follows 

S = J d 4 x^^^+^(^ 1 (a)i? + G^ 2 (a)^ + C^J- 4 (n)C^) - a) (4.2) 

where we have used that Gy is traceless and have defined J r 1 (n)+ij r 2 (D) = .Fi(D). Equations 
of motion for action (4.2) can be derived by substituting Ry = Gy + \5yR in (2.2) 

MpG* = T£ - A<5£ - 2\G£f 1 (n)R + 2A(V^ - 6»n)fi(n)R - l\Rr 2 (p)Gr- 



-2XG^ 2 {U)& V + ±6ZG%JF 2 (u)& a + 

+2A (v p V v F 2 {n)&r - ^OF 2 (n)G£ - a V p F 2 (n)G^ + 

- ^ {££ + A) + A£ 2 £ - ±% {££ + £ 2 ) + 2AA£ + 2AC£ 
where we have defined: 

oo n—1 oo n—1 



(4.3) 



(n-0 

Ji -it- 

n=l Z=0 n=l i=0 

n—1 oo n—1 



W = E /2n £ V^)«V„G^-^, £ 2 = E ^n E G®^. 

71=1 «=0 71=1 i=0 

oo n—1 

A£ = £ hn E V^G^V^G^" 1 ^ - V^G^G^-'- 1 ^] 

71=1 £=0 

and fi n are coefficients of the Taylor expansion of function T\. The Weyl tensor related part 
may have an impact now since we are going to consider perturbations. One can find the 
relevant part of it is 

C£ = (R aP + 2V Q V /3 ) T 4 (n)C u a ^. 

Saying relevant we mean only the piece which is obtained by the variation of one of the Weyl 
tensor factors in the action. Then we are left with only one Weyl tensor as it is obvious from 
the latter formula and further perturbation of this remaining Wei tensor factor may produce 
a non-zero contribution to the perturbation equations. 
The trace equation becomes 

-MpR = T - 4A - 6\nfi(\3)R - A(A + 2£~i)- 
-2W p V n F 2 {u)& lp - A(£ 2 + 2£ 2 ) + 2AA 

and the Weyl tensor related term C turns out to be zero thanks to the full tracelessness of 
the Weyl tensor. 

This form of action and equations of motion also turns out to be beneficial for studying 
perturbations as we will see in the next Section. 
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4.3 Avoiding recursion relations 

Surely we are still sticking to the FLRW type of the metric but keep up to some extent the 
scale factor general. This means that we cannot say anything specific about functions 
which appear after the box operator acts on the scalar curvature n times. Fortunately, we 
now can have a tiny progress with the second rank tensor Gu- It is traceless and action of 
the box operator does not break this property because box commutes with the metric. It 
means that = for any n and thus we can introduce 

U n &t = b n {t)C where # = diag(3, -1, -1, -1) (4.5) 

where one can compute that 

b n+1 (t) = (□ + 8H 2 ) b n (t) and b (t) = H/2. (4.6) 

Here as usual H = a/ a is the Hubble function and dot is the derivative with respect to the 
cosmic time t. 

This allows us to simplify all the tensor structures in the trace equation (4.4). As was 
explained in the previous Section this is the cornerstone equation since a solution to it is 
almost automatically a solution to all the Einstein equations modulo perhaps some radiation 
(assuming no other matter is present in the system). Careful substitution gives 

oo n— 1 

C 2 = 12 hn (# M ¥Wn-i-l(i) " aff^&n-t-lW) , 
n=l 1=0 
oo n— 1 

£ 2 = 12^/ 2n ^(*)&n-z(i), 

n=l 1=0 
oo n— 1 

A = -48^/ 2n ^ (ff0 t (6/(*)&n-z-i(t)) + (3# 2 + #)M*)&n-;-i(i))) , 

n=l 1=0 

oo 

V p V„T 2 (a)G^ = 3ab(t) - 12H(b(t) + Hb{t)) where b(t) = ^(/ 2n 6 n ). 

71=1 

(4.7) 

It is manifest now that we have reduced the problem of solving the trace equation to the 
problem of many scalar functions and moreover all of them are connected by means of the 
covariant box operator. To have the complete picture we note down 

R = 12H 2 + 6H. (4.8) 

Possible recursion relations may be of course substituted here. Technically presence of a 
recursion relation for w- n ^ or b n or both means that the corresponding infinite series recurses 
after a certain n. The following configurations may appear than in this model: 

1. If both and b n obey some recursion conditions for a given a(t) then this in turn 
means that there is a finite set of linearly independent functions and in order to satisfy 
the trace Einstein equation we must equate to zero coefficients in front of each such 
function. For a finite n this produces finite number of conditions on /i, 2(t - This is 
exactly what is implemented in the previous Section for the solution (3.2). 
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2. If only, say, obeys a recursion relation while b n not then this means that in general 
infinite fine tuning is required. Indeed, absence of recursion for b n is equivalent to the 
fact that infinitely many linear independent functions will be generated all accompanied 
by fi n coefficients. Here we see two ways to satisfy the trace Einstein equation. Either 
all the coefficients f2 n must be trivial eliminating the corresponding action modification 
at all or we must fine tune all of them such that the infinite series of functions sums up 
to a constant. The latter possibility may be accomplished in the spirit of the so called 
"addition theorem" for the associate Legendre functions when some certain infinite 
series sums up to the unit (see for instance formulae (8.794) and (8.814) in [39]). For 
this to happen it is vital to have infinite number of non-zero coefficients fi n . This is 
reasonable in our set-up but on the other hand all of them must be fine tuned. This is 
not definitely unacceptable, however. When, for example, we speak about a function 
like J~2 = eP u we say it has one free parameter j3 rather than infinitely many. The 
point here is that these operator functions come normally from the consideration of 
the corresponding string models as mentioned in the Introduction and their nature is 
therefore originates from some other principles and not from our cosmological model. 
Thus in the case some solution needs an infinite number of coefficients to be fixed we 
should only worry whether the resulting operator function is acceptable from the point 
of view of the string model it comes from. Example with such properties is solution 
(3.3). 

3. If neither nor b n appear to obey a recursion all the coefficients /i,2 n are in the 
play. Again they can be either trivialized or fine tuned. Here, however another type 
of fine tuning may arise when instead of fixing all the coefficients we do fix only some 
relations in between of f\-s and /2-s. This is still an infinite fine tuning from the point 
of view of individual coefficients but from the point of view of the operator functions 
J- this may be just one simple relation. To be more precise, provided cancellation of 
all terms happens if /i n //2 n = const which does not depend on n then it is equivalent 
to have T\ ~ Ti. 

4.4 Other solutions 

Solutions of the first type. 

Solution belonging to the first type when both the scalar curvature and the Ricci tensor 
obey recursion relations, apart from (3.2), is 

a(t) = a t (4.9) 

One can find that on this solution 

R = 6/t 2 , OR = 0, 6 = , h = -~, 6 2 = 

and it is easy to show that equations of motion can be satisfied. It is not a bouncing solution 
but may be interesting in analyzing the regime when the scale factor grows linearly. 

Configurations of the second type. 

Known configurations of the second type are as follows 

a(t) = a Vi (4.10) 
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For this scale factor we have 

113 

R = 0, b = h = —, b 2 = . . . 

and sequence of b-s does not stop. 

Equations become a bit more transparent if we pass to a(t) as the variable. We can do 
this in each region where we have one to one correspondence t a(t). We therefore can 
rewrite equations in the new variable and time derivatives transform to dt = aH(a)d a . 
The most important for us d'Alembertian operator □ becomes 



a H 

where prime is the derivative with respect to a. Another important operator is 

□ +8 ^=-^(^ + (i + |;) 8 .-A 

since it creates functions b n . 

To have a connection with the previously obtained solution note that 



1-^| 



gives a(t) = aocosh(at). 

Than one can guess some interesting configurations and one of them is 



1-4 



giving 



a(t) = ao^/cosh(at) (4.11) 

We have for this scale factor 

R = 3a 2 

while 6-s form an infinite sequence of polynomials of cosh(cri). 

One interesting observations is that the scale factor given by (3.3) is exactly the leading 
term in the series expansion 

H = a\f^ 

V dP 

in the limit p — > 0. Also note that for scale factor (3.3) we have 

R = I2a 2 t 2 + 6a, DR = -72a 3 t 2 - 24a, . . . 
b = a/2, h = Aa 2 t 2 , . . . 

where we see explicitly that R obeys a recursion relation while b not. 

In all these cases one need an infinite fine tuning as explained in the previous Subsection 
in order to perform summations in the equations of motion. Even though such a fine 
tuning may be not favorable the reformulation itself of equations of motion in terms of 
a may be interesting and we are looking forward to go further in this direction. 
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Note also that all the above configurations of the second type are solutions with Ti = 0. 



Configurations of the third type. 

The third type of configurations is extremely generic and the only known scale factor 
which makes us able to track equation up to some reasonable extent is 

a = a t p (4.12) 

for a generic p. For this scale factor we have 



n 



' 1 j-21+2 



where constant coefficients may be expressed through p explicitly. Note that for odd p the 
series for R l stops at a certain I and the soluions falls back in the previous subclass when one 
recursion relation exists. Now it is a matter of a straightforward computation to see that 
one can satisfy the Einstein equations by imposing algebraic relations on coefficients f\ n and 
/2 n . This is because internal summations over I in expressions for C-s do not involve the 
time variable t. 



Note also that all the above solutions are still solutions with an arbitrary T±. Even 
though it is not important for the background one would change the perturbative picture 
playing with the JF4 parameter. 



5 A road- map to perturbations 

Needless to say that even writing down the perturbation equations may be an unfeasible task 
in the above formulation of the model. Below in this Section we outline the general approach 
and emphasize what shall be developed shortly in subsequent papers since the full detailed 
analysis is expected to be rather cumbersome and deserves a separate paper. 

5.1 Bianchi identity 

We know that the Einstein equations are constrained with the Bianchi identity which says 
VuG{! = 0. In our case we have more than just Einstein-Hilbert Lagrangian but all the 
additional ingredients we have are covariant terms. This guarantees that the Bianchi identity 
holds trivially without imposing any extra condition. On the other hand this implies thanks 
to arbitrariness of coefficients fi n that each separate term does covariantly conserve. Indeed, 
each fj n is a coefficient in front of some covariant structure in the Einstein equations, say 
t„ . Assuming that only one of f-s coefficients is non-zero we come to a conclusion that the 
corresponding structure must covariantly conserve due to Bianchi identities, i.e. V^r^ = 
0. In other words it resembles a conserving perfect fluid stress-energy tensor. The same 
argument is applicable to all the f-s coefficients as well as their arbitrary combinations. 
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The above arguments imply that thanks to Bianchi identity the parts with different Ti 
covariantly conserve separately To make use of this we define 

Tot = T£, 

T\t = -2AG{J Ji(D)i2 + 2X(Wd u - ^D) Ji(D)i2 + XC^ - (££ + A) , 



+ 2A ( V p V u F 2 {u)G»p - tuT 2 (U)Gt - CT V P T 2 {U)G^ ) + 



l -\RF 2 {u)& v -2\G»F 2 {u)& v + K 
l -UT 2 (u)GZ- l -< 

+ A£ 2 £-^(£ 2 £ + £ 2 )+2AA£, 
T 4 £ = 2A (R a p + 2V a V /3 ) F±{n)C u a ^. 
Now Einstein equations can be written in an extremely concise form 

m^ = ]Tt^-a^ (51) 

and moreover we have 

V^Trf = for any /. 

One recognizes here the system of minimally coupled perfect fluids minimally coupled 
to gravity. The perturbation technique is known in general (see [32, 40] for instance) but it 
is not obvious it can be applicable "as is" to our model as we will see shortly. 

5.2 Subtleties of the non- localities 

Having brought Einstein equations in our model to an already studied structure does not 
mean we are able easily, if at all, solve the perturbation equations. 

It looks promising that we could clearly separate the whole equations into minimally 
coupled terms. But we have to stress here that each term is not a canonical perfect fluid, 
i.e. it cannot be written just in terms of energy density p, pressure density p and four-fluid 
velocity n M as discussed in [41]. 6 This results in presence of anisotropic stresses as well as 
entropic perturbations for our fluids at the level of perturbations. This in turn does not allow 
us to use the already known techniques straightforwardly since the system of perturbation 
equations (even non-local ones) will not be closed. 

On the other hand all Tj are not independent external stress energy tensors but rather 
some structures built on the metric. We therefore must be able in principle to have the 
corresponding perturbative quantities in terms of the metric perturbations. 

Also one must not go straight with T4 term since it represents a contribution which is 
absent in the background. It means, that energy and pressure for T4 are both zero. Thus it 
is not obvious how one can define the equation of state parameter w = p/p and the speed of 
sound c s 2 = p/p. 

It is, however, the simplest term for performing actual calculations at the perturbative 
level. This is because it is zero on the background and completely traceless. The first point 
means that 



5T 4 » = 2A {R a p + 2V a V /3 ) F A (U)8C V 



a/3fi 



3 In other words we cannot just have Tiy = (p + p)u M ti„ — p&t- 
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The second point means that for isotropic perturbations 

where rj is the conformal time, k is the comoving wave- vector and P u Q ^ is a constant tensor 
which keeps all the symmetry properties of the Weyl tensor. We thus can use the same 
idea as in the previous Section for the traceless tensor G„ and compute an action of some 
differential operator T> on 8C V as 

V5C v a ^ = {Vc(j])y U 'P u a/3fM 

with a simple algebraic relation between T> and T>. 

We can push these arguments further and draw some conclusions when R or Gy are 
perturbed in T\ and T 2 . We then are able to compute implicitly action of the non-local 
operators on those perturbed quantities. However, the real problem comes would one perturb 
the non-local operators themselves and compute the action of this perturbed operators on 
unperturbed background quantities. This is where presence of recursion relations is crucial. 
This was the key point why could we end up with final results in [18]. It is an open question 
how one would go in a general case when there are no recursion relations. We hope to address 
this issue in the forthcoming research [42]. 

We nevertheless can come up with a conclusion that the discussed above solution (3.2) 
is an example of a stable bounce configuration. Even though this particular solution creates 
recursion relations and in principle one would much easier write down system of perturbation 
equations we can claim the bounce phase is stable without tedious calculations. This stability 
is guaranteed because solution has analytic dependence on coefficients /. Stability of this 
solution with T 2 = JF& = was explicitly demonstrated in [18]. Therefore turning on JF 2 
and/or J-4 must keep the configuration stable up to some range of coefficients /2,4 n - Of 
course, this argument is not enough to find the allowed domain for new coefficients and this 
is the primary goal for the forthcoming paper [42]. 

5.3 de Sitter limit 

The de Sitter limit is perhaps the most simple configuration for analysis of perturbations 
(apart from Minkowskian space-time, of course). The de Sitter Universe is described by the 
scale factor 

a = a e Ht 

and the corresponding Hubble function is just the constant H. This is where we see why 
the traceless tensor G> is convenient since it is identically zero in such a metric. Moreover 
the scalar curvature is a constant R = 12H 2 . One can easily check this is a solution to 
our equations. But the most important that the part of equations which may be relevant 
for perturbations is much shorter since we can drop terms quadratic in Gy as well as some 
others. Then the relevant equations become 

M 2 P 5G^ = ST? - 2XSG^f l0 R + 2A(V^ - ^□)Ji(D)J J R - hjUF 2 {U)8G$- 

+2A \V p V v F 2 (p)6Gi* - ^UF 2 {U)80t - a V p P 2 (a)8G^ - (5.2) 

~8£R (ji(D) - / 10 ) SR + 4XV a V^ 4 (a)5C u a ^ 
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with the trace 

-M%8R = 6T- 6ADJi(D)5i2 - 2XR (ji(n) - / l0 ) 5R - 2AV p V M J" 2 (n)<5(5 w . (5.3) 

These can be analyzed up to the end explicitly for some configurations. To highlight the 
track towards solutions to the perturbation equations for scalar perturbations we note that 
5R has an explicit expression in terms of Bardeen potentials <3? and Vf, the two gauge invariant 
scalar degrees of freedom [43]. Then one can show that i ^ j component of equation (5.2) is 
a non-local equation containing 5R and <J> — ^, where i,j are the spatial indexes. Assuming 
further that the matter is the radiation without anisotropic stresses and computing explicitly 
VpV At J r 2(n)(5G w which is feasible thanks to the tracelessness of the G^ p tensor we must end 
up with two non-local equations on <I> and . 

Even though we can develop on this right away we are going to stop at this stage because 
there are many other serious and unrevealed questions regarding the de Sitter background 
in this type of models such as the ghost-free condition. These unexplored matters deserve a 
separate study which will include the detailed analysis of perturbations as well [42] . 

6 Summary and Outlook 

We have considered in this paper the most general extension of GR based on inclusion of 
stringy motivated non-local operators and keeping the quadratic in curvature terms. 

The primary goal was to find out an exact bouncing solution and it is for the first time as 
this task is accomplished up to the best of our knowledge. It is intriguing that the bouncing 
configuration is the cosine hyperbolic which is still a solution provided operator functions 
J~2 and J-4 which control the presence of extra terms are trivial. This indicates a presence 
of some symmetry one would have to find and so far is one of the open questions. A good 
point of our bouncing solution (3.2) is that it may be a solution in our model without extra 
matter which was not the case without Ti and J-±. Also this solution requires just a couple 
of conditions which are very general and are easy to satisfy. 

Moreover we have presented several other solutions which are not all bouncing solutions 
but anyway widen our understanding of the model. Practically we divided possible solutions 
into three groups based on the fact whether they create or not recursion relations when 
the scalar curvature and the Ricci tensor are being acted by powers of the d'Alembertian 
operator. Among all the found configurations the one which is (4.11) looks very interesting 
for the further study since it creates a bouncing Universe with the constant scalar curvature. 
It is a solution with T% = or for a very special Ti when all the coefficients are strictly fixed. 

One more interesting point of this model is that the last term with the Weyl tensor 
squared in the Lagrangian does not contribute to the background at all. It, however, shows 
up in perturbations and it is interesting to see whether presence of this term can be efficiently 
used to tackle the problem of solving the perturbation equations. 

Also we were able to bring the model formulation to such a form that was already 
analyzed at the perturbative level for other configurations. The perturbations can be analyzed 
in full at least numerically in the de Sitter limit and the clear way towards this was outlined. 
As it was mentioned in the main part of the paper, however, we postpone the more detailed 
and technical study of this question since it is more natural to join it with other unrevealed 
problems one can put for this model in the de Sitter space. Namely, ghost-free condition and 
subsequent quantization are the issues to be fully studied as well and results should appear 
shortly in [42]. 
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The full analysis of perturbations is expected to be very hard in general. At the moment 
it is not even obvious one will be able to close the system of perturbation equations for 
scalar perturbations. Tensor perturbations may become simpler but still one would face the 
difficulty that effective anisotropic stresses present in the system. This is tough and open 
question at the moment. 

Furthermore one question which we kept aside is the anisotropic perturbations during 
the contraction phase since this question is under investigation in the parallel project [44]. 
The generic expectation is that such perturbations must grow during the contraction phase 
and the main problem is to formulate the domains of parameters of the model allowing the 
bounce to happen smoothly. Approaches to study anisotropic perturbations in non-local as 
well as stringy-inspired models can be found in [22, 45] and references therein. 
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